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Abstract 



We discuss some 3D numerical simulations of accretion discs using the 
SPH method and a polytropic equation of state. We show that discs exist even 
for as large value of the polytropic index as 1.2, and that these discs are always 
in hydrostatic balance. We also show that even without any inflow, spiral shocks 
Q-f appear in the discs. 
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+3 . 1. Introduction 

CO 

Numerical studies of accretion discs have been mostly restricted to 2D 
cases, due to computing time limitations. Among many things, these 2D sim- 



ulations have shown that spiral shocks appear in inviscid discs (e.g. Sawada et 
al. 1987). Recently some 3D simulations have been carried out (see Yukawa, 
Boffin & Matsuda, 1997 for an uncomplete list), mostly using particles methods. 
These simulations were apparently unable to generate spiral shocks in the accre- 
tion disc, but this could be related to the fact that they used either an isothermal 
or pseudo-isothermal equation of state, either neglected pressure effects or used 
too low resolution. We have run three-dimensional Smoothed Particle Hydrody- 
namics (SPH; see e.g. Monaghan 1992 for a review) simulations with a polytropic 
equation of state. This method includes self-consistently the effect of pressure 
forces and we checked that we could always resolve the disc in the vertical dimen- 
sion. Concerning this last point, we therefore used a variable smoothing length 
(which, in SPH, decides the resolution) and checked that at each point in space, 

* To be published in Disk Instabilities in Close Binary Systems - 25 Years of the Disk- 
Instability Model, Mineshige S. & Wheeler J.C., eds, Universal Academy Press, 1999 
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the smoothing length, h, was smaller than the disc scale height, H. Details of the 
method and of some of the results can be found in Yukawa et al. (1997). 

3D; 7 = 1.2; 59000 ptcls 
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Fig. 1. Structure of the flow in the orbital plane, for our simulation with a 
mass ratio of 1.0 and 7 = 1.2, after 2 orbital periods. Only the region close to 
the primary is shown. A total of 59,000 particles was used. 

2. Results : mass inflow 



In figure |T], we show the flow at the end (i.e. two orbital periods) of our 
simulation with mass inflow when we use a polytropic index, 7=1.2 . As can 
be seen, a spiral structure is clearly present, confirming the fact that SPH is 
able to tracks these structures but, more importantly, that these structures are 
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present in 3D accretion flows. This result also confirms that a disc does form in 
3D, even for such a large value of the polytropic index. Moreover, the disc is in 
hydrostatic balance, as its disc height is precisely equal to the value expected: 
H ~ c s /Q, where c s is the sound speed and fl is the angular velocity. Because, we 
use a rather large sound speed as initial condition (0.1, where the orbital velocity 
corresponds to 1.0) and a large polytropic index, the disc we obtain is rather hot, 
hence rather thick (— ~ 0.2 — 0.3). For the smaller value of 7 used, 1.1 and 1.01, 
we obtain smaller disc heights : 0.12 to 0.2 and 0.09, respectively. In both cases, 
the hydrostatic balance in the vertical direction holds true. And in all cases, 
the ratio between the local vertical disc height (i.e. the disc semi-thickness) and 
the local smoothing length lies between about 2 and 6. Thus, we have certainly 
resolved the disc vertically. 

Just a note in passing concerning the viscosity present in our code. We use 
the standard artificial viscosity of SPH which, as shown e.g. by Murray (1996), has 
an equivalent shear viscosity, v ~ 0.1 aspnc s h. In term of the Shakura-Sunyaev 
a- viscosity, ass, this can be rewritten, 

«sph ^ h a S PH h 

ass — — • 1 

SS 10 c s 10 if K 1 

With the value of «sph = 1 used, we therefore have an equivalent ass of 0.02 to 
0.05. 

It has to be noted that we cannot claim to have obtained a true steady 
state as the mass in the disc is still increasing at the end of the simulations. 
Two-dimensional simulations (Boffin et al., in preparation) show us that several 
tens of orbital periods are necesary to reach a steady state. However, in our 3D 
simulations, we can see that the structure of the flow does not change after, say, 
one orbital period. We therefore believe that we have reached a "quasi-steady 
state" and can study the final structure of the flow. We cannot, however, make 
any definite claims about the mass accretion rate. From figure [l], we also observe 
that we do not have a true "hot spot" but more a kind of "hot line". This is, we 
believe, again due to the large initial sound speed, resulting in a very wide inner 
Lagrangian stream. 

In figure 0, we show the same as in figure [L], except that we have divided 
the particles following their height above the orbital plane. This can be used to 
study the possible variation of the disc height with the orbital phase as obtained by 
Hirose et al. (1991). We do not seem to find any conclusive variations, however. 
Also, we cannot observe any stream overflow in the z-direction as obtained by 
Armitage & Livio (1996). The reason for this discrepancy is unclear and we are 
presently working on this. Possible reasons are : their use of a large viscosity, 
their initial conditions, our large initial sound speed, ... 
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3D; 7 = 1-2; 59000 ptcls; sliced in 13200 ptcls 



< |z| < .0075 


.0075 < |z| < .0175 


.0175 < z < .0335 


.0335 < |z| < .070 



Fig. 2. Same as fig.|T], except that we have subdivided the particles following 
their height, z, above the orbital plane, in such a way that there are as many 
particles (13,200) in each frame : 

a) 0.0000 < z < 0.0075, 

b) 0.0075 < z < 0.0175, 

c) 0.0175 < z < 0.0335, 

d) 0.0335 < z < 0.0700. 
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3D; 7 = 1.2 and 1.01 ; q=0.5; 9,000 ptcls; T=1.39 




3D; 7 = 1.2 and 1.01 ; q=0.5; 30,000 ptcls; T=2.00 



Fig. 3. Results of our simulations without mass inflow for two value of the 
polytropic index, 7=1.2 (left) and 1.01 (right). The upper frames correspond 
to a lower resolution (9000 particles) than the lower frames (30 000 particles). 
The mass ratio is 0.5 and the plots represent greyscale maps of the asymmetric 
part of the density field. Note the effect of increasing the resolution. 
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3. Results without mass inflow 

We have also performed several simulations without any mass inflow. In 
this disc is initially set-up around the primary, so that it is resolved 

vertically and in hydrostatic balance. It is then evolved with the full potential of 
the binary system taken into account. Here again, as shown in figure [J], which 
is a greyscale map of the asymmetric component of the density, spiral shocks can 
clearly be seen, both in the 7=1.2 and 7=1.01 cases. Thus, these spiral shocks are 
not the result of the inner Lagrangian flow. This is not a surprise if, as believed, 
the spiral structures are due to the tidal force of the companion (e.g. Savonije et 
al. 1994). 

Figure also shows the importance of resolution : although with 9,000 
particles we cannot find any severe difference between 7=1.2 and 1.01, this is 
no more true with 30,000 particles. For 7=1.01 indeed, in the inner part of the 
disc, the spirals become more tightly wound, a result well known in 2D (e.g. 
Sawada et al. 1987). The reason for this difference may lie in the fact that for the 
7=1.2 case, the Mach number of the flow always remains smaller than 10, while 
for the 7=1.01 case, it starts at a little below 10 in the outer part of the disc 
to reach above 30 in the inner part. It was already shown by, e.g., Savonije et 
al. (1994) that the higher the Mach number, the more tightly wound the spiral. 
What is noticeable in our 3D simulations, is the fact that we cannot make any 
clear distinction between a cooler disc (7=1.01) and a hotter disc (7=1.2) when 
we restrict ourselves to the outer parts of the disc. If confirmed, this result may 
be useful to reproduce the observations of spiral shocks in IP Peg (Steeghs et al. 
1997) even when we deal with the cool disc typical of a dwarf nova. 
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